Quantum simulation of an extra dimension 
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We present a general strategy to simulate a D+ 1 -dimensional quantum system using a £)-dimensional one. 
We analyze in detail a feasible implementation of our scheme using optical lattice technology. The simplest non- 
trivial realization of a fourth dimension corresponds to the creation of a bivolume geometry. We also propose 
single- and many-particle experimental signatures to detect the effects of the extra dimension. 



Introduction. There is long-standing interest in Physics for 
the possible existence and effects of extra dimensions. This 
interest was brought by the seminal papers by Kaluza and 
Klein [1] aimed at unifying interactions using the presence 
of a fourth spacial dimension. Later on, extra dimensions be- 
came a sine qua non element for the construction of string 
theory [2]. Separately, enormous progress has been made in 
recent years to achieve real quantum simulations, that is, to 
simulate quantum mechanical models using other well con- 
trolled quantum systems [3]. It is now reasonable to investi- 
gate to what extent present technology can be used to faith- 
fully simulate a quantum theory living in extra dimensions. 

Let us briefly recall recent progress on quantum simulation 
of condensed matter models using cold atoms [4, 5], By con- 
fining atoms to an optical lattice, the Hubbard model may 
be realized [6], and the superfluid-to-Mott-insulator transi- 
tion observed [7]. Furthermore, several schemes to couple 
neutral cold atoms to artificial Abelian [8] and non-Abelian 
[9, 10] magnetic and electric fields have been put forth [11- 
14]. This opens the door to creating strongly correlated 
quantum-Hall states with cold atoms [15]. Although cold 
atoms are non-relativistic, it is possible to simulate relativistic 
effects by looking at the low-energy behaviour of some spe- 
cial lattice models -lattice models where the band structure 
presents Dirac cones, e.g. honeycomb lattices [16, 17] or lat- 
tices dressed with internal degrees of freedom [18]. Hence, it 
is not far-fetched that that cold atoms may provide some in- 
sight into particle physics models that are not completely un- 
derstood, such as quantum chromodynamics (see for instance 
[19]) or, as presented here, in the analysis of extra dimensions. 

General strategy. The basic idea to achieve a quantum sim- 
ulation of an extra dimension consists in engineering the con- 
nectivity of the system partly on real dimensions, and partly 
on the use of different species for the degrees of freedom 
(for an alternative approach cf. [20]). Let us illustrate this 
construction in the case of the simplest quantum mechanical 
model of a free particle on a hypercubic D+l spacial lattice. 
The Hamiltonian for this system is 
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ity of the lattice. Now, we write the D + 1 -position index q as 
the combination of a D-dimensional position r and a separate 
index a in the extra dimension. That is, the D + 1 -space is de- 
composed in hypersurfaces, that we shall call layers, labeled 
by the a index, q = (r, a), and the Hamiltonian becomes 
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H.c. 



(2) 



where a q destroys a degree of freedom at site q, and the u 7 
stand for D+l cartesian versors that set the actual connectiv- 



The operator a 1 - ' can be reinterpreted as a Fock operator 
for the species a. The total number of species N, <J = 
l,...,N, corresponds to the number of lattice layers in the 
extra-dimension. The crucial requirement to simulate a D + 1 - 
dimensional model with a set of N species in a D-dimensional 
lattice is that each internal state be coupled to only two other 
states in a sequential way. 

We here propose two ways of realizing the Hamiltonian of 
Eq. 2. Both methods make use of Raman transitions in optical 
lattices but differ in the way of constructing internal degrees 
of freedom. We shall refer to the two methods as: 
i) State-dependent lattice; if) On-site dressed lattice. 

State-dependent lattice. Let us now discuss the realization 
of our model in an extra dimension using spin-dependent lat- 
tices [21] (alternative state-dependent lattices could be ob- 
tained by trapping long-lived optical states via superlattice 
techniques [17, 22]). The simplest non-trivial step to construct 
an extra dimension corresponds to a bivolume geometry, that 
is, the quantum system spans over N = 2, D ~ 3, layers con- 
nected through an extra dimension. We shall keep this exam- 
ple in mind, while providing expressions for N generic layers 
in D dimensions. 

It is possible to induce a relative phase in the periodic opti- 
cal potential seen by hyperfine states with different total angu- 
lar momentum F, by tuning the angle between the linear po- 
larizations of counterpropagating lasers that form the lattice. 
Hence, the respective minima seen by the atoms get separated 
by the same phase shift. For the explicit example of &7 Rb 
loaded in a 3D lattice, the two available values of F, F = 1 
and F = 2, become two different species, which correspond 
to the two layers of a bivolume. In general, if we have N dif- 
ferent minima and the atomic states are long-lived, the system 
will consist of N copies of 3D layers. 

At this stage, there is no hopping between different layers. 
The additional ingredient required to map uncoupled 3D hop- 
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Figure 1. A bivolume geometry is made out of two 3-dimensional 
sublattices separated by a displacement vector d. Lattice sites in dif- 
ferent colors trap different internal states. Black and red links con- 
nect nearest neighbors of 3D sublattices, respectively, and represent 
the free hopping terms in each sublattice. These tunnelings are due to 
kinetic energy and do not involve transitions between internal states. 
Green links connect nearest neighbor sites from different sublattices 
and are induced by laser assisted Raman transitions. The strength of 
this coupling depends very strongly on the distance between the pairs 
of sites. 

ping Hamiltonians to a single AD Hamiltonian is a coupling 
between different minima or layers. The main difficulty is to 
ensure that each lattice site have the appropriate number of 
nearest-neighbours. In a regular optical lattice this is always 
the case as Wannier functions -single-particle functions liv- 
ing on each lattice site- at different positions are orthogonal. 
To analyze this requirement, let us order arbitrarily the hy- 
perfine states and label them with the index a = 1 . . . , N. We 
now consider a Lambda configuration between the state a and 
the state a + 1 which induces a Raman transfer between the 
two species, and hence an assisted hopping term between the 
two sublattices. The effect of the laser is shown in Fig. 1 for 
D = 3, where for simplicity only one cubic cell of the two sub- 
lattices is pictured and the spectator species are omitted. The 
two sublattices are separated by a displacement vector d, de- 
fined as the smallest distance between a pair of vertices of the 
two cubes. Black and red links represent free hopping terms 
within the sublattices while green links correspond to Raman 
transitions. In order to make the picture readable only the 
stimulated hoppings relative to only one site in the a + 1 sub- 
lattice are displayed. For each of the possible hopping terms, 
the Raman hopping rate is computed as the overlap integral 

Jal+i = ^Jd D r w*(r - />(r - d) , (3) 

where the V , j = 1,... ,2D, are the positions of the blue ver- 
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Figure 2. Behavior in log scale of the hopping hierarchy under 
changing of the maximal Cartesian component d x of the displace- 
ment d for different values of the lattice potential V. From top to 
bottom, the curves correspond are to V = 5, 10, 15, and 20Er. The 
dash horizontal lines give the suppression by a factor 10 and 100 of 
the next-to-leading order coupling. Small values of R allow for the 
proper bi-volume connectivity of the system. 

tices measured with respect to the front-left-down one, 
and w(r) is the Wannier function centered at r, which for a 
hypercubic lattice is the product of the 1 -dimensional Wan- 
nier w of each Cartesian direction. In a more complicated 
scenario, as superlattices, the Wannier functions for the two 
species may be different. 

Wannier functions are localized and decay exponentially 
fast away from each lattice site. A hierarchy between the hop- 
ping rate /^ct+i an< ^ tne -^al+v J > ^, can be easily generated, 
already for not so small value of \d\/a, with a the lattice spac- 
ing. The behavior of the suppression depends on the depth of 
the optical potential V, and increases with V. By using the 
separability of the Wannier functions, the maximal ratio be- 
tween the unwanted links and the first link can be computed 
as Id problem 

= Max{j^ +lJ >2} = Jdxw*(x)w(x + a-d x ) 

JdxwXx)w(x-d x ) ' ; 

where d x is the largest Cartesian component of displacement 
d. The optimal scenario at fixed |d| is for d along the diagonal 
of the hypercube, i.e. d x = The ratio R is plotted in log 
scale for typical values of V about ten times the recoil energy, 
Er = p^2) . in Fig. 2. The suppression is very efficient: for 

V = 20Er and d x = | the other hoppings are less 1 % of /^i+l ■ 
Under this condition, the effective Hamiltonian is 

H = -l^{tja[%a^+J'af +X) "af ) ) + H.c., (5) 
r,a ;'=1 

where for simplicity we assume a uniform hopping between 

species J' = J^ +l , Vcr. The above Hamiltonian coincides 
with the one of Eq. 2, and hence is equivalent to the free 
hopping Hamiltonian Eq. 1 in D + 1 -dimensions. 
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A few considerations are in order. First, we may choose 
periodic or open boundary conditions for the hopping term 
in the extra-dimension, by either including or not a Raman 
stimulated transfer between the N and 1 states. In the for- 
mer, our model is equivalent to D + 1 -model compactified 
on a circle (cf. [1, 2, 23]). Furthermore, as non-zero con- 
stant or r-dependent phase can be chosen for the hopping rate 

the above set-up allows the simulation of compatifications 
in presence of non trivial back-ground magnetic fluxes pierc- 
ing the circle. Second, if the interactions are negligible and 
open boundary conditions are chosen, the only limitation to 
the "tickness" of the extra-dimension is due to the number 
of hyperfine states available, and to the technical difficulty of 
coupling them selectively. Indeed, for any value of potential 
V we can find a displacement d sufficiently small such that 
R <C 1 and Eq. 5 holds. 

Once quartic interactions due to binary scattering of atoms 
are included, the non-zero overlap of the Wanniers of different 
D-dimensional sublattices potentially induces nearest neigh- 
bours interactions in the D+ 1 -dimension. Again, the sepa- 
rability of the hypercubic lattice allows to express the ratio 
between the nearest neighbours and onsite interactions as 

2 

i B s _ <Wn -fi Jtej \w(xj)w(xj-dj)\ 
(Ru) aa+ i = — [I . |4 ' W 

"CTCT j = \ J dXj W(Xj) 

where dj, j = l,...,D, are the Cartesian components of the 
displacement d, and a aa+ \ and a a a are the scattering lengths 
of <7-(T+ 1, and <7-<7 collisions, respectively. In most cases, 
Woo+i ~ cc aa , hence (Ru)aa+i = Ru- in the Table I the ratio 
Ru for a displacement along the diagonal of the hypercube, 
d x = dj, \/dj, is given for different values of V. As expected, 
Ru is monotonically decreasing function of both d x and V. 
By comparison of the numerical result of Table I with the plot 
2, it turns out that nearest neighbours interactions cannot be 
disregarded for small value of the optical potential (V = 5Er) 
while the required hierarchy for both hopping and interaction 
couplings is realized for a range of displacements d x at larger 
values of V. In fact, a non-trivial Ru may lead to interesting 
new phenomena like supersolidity [24] . 

On-site dressed lattice. In another way to realize of Eq. 2, 
the internal atomic degrees of freedom can be obtained con- 
sidering hyperfine states with same total angular momentum 
F and different third component nif . In our case, we are in- 
terested in atoms with sufficiently large values of F, in order 
to have several layers in the extra dimension. Very recently, 
earth-alkali atoms with such property have attracted a lot of at- 
tention from a theoretical [25, 26] and an experimental point 
of view [27-29]. As in the ground state they have total orbital 
momentum J = 0, and the nuclear spin / practically does not 
couple to the dynamics, it follows that F =1, and their interac- 
tion are S\J(N) invariant with N up to 10 for 87 5V. Such atoms 
are fermions. Smaller symmetry groups can be achieved also 
for alkali atoms [30]. The mixing interaction between species 
can be realized by optical means as in [31]. It is worth notic- 
ing that there is no additional difficulty in: i) realizing periodic 



boundary conditions by identifying the species N+l and 1 , as 
this amounts to coupling species 1 and N as well; it) engineer- 
ing hoppings with non-trivial phases equivalent to magnetic 
fluxes piercing the compactified circle. If 2-body interactions 
between species are present, the corresponding extended Hub- 
bard model has a very rich phase diagram that is at present 
under study [32]. Such scenarios are technically very chal- 
lenging and involved but potentially feasible using state of the 
art techniques. 

Single-particle observables. We now turn to discuss ob- 
servables that can reveal the properties of an extra dimension, 
whatever experimental technique is used to implement it. The 
most direct way is to consider single-particle effects. The ef- 
fective dimensionality of a system is revealed by the scaling 
behavior of observables. For instance, the two-point correlator 
free bosonic relativistic field theory in D spacial dimensions, 
decays as e -m W/|jc|^ _1 \ where \x\ is a space-time distance, 
m is the mass of the field. Such dimensionality dependence 
can be interpreted as the effect of monogamy of entanglement 
[33]. The more neighbors exist, the more distributed the cor- 
relations must be. In this simple framework, the effect of an 
extra dimension which can be adiabatically coupled in the sys- 
tem must translate into an interpolation between the decay ex- 
ponents in the propagator. This is managed by a tower of the 
so-called Kaluza-Klein modes that bring a series of massive 
states into the spectrum. We show here that the local density 
of states, p(E), routinely measured in cold atoms experiments 
[4, 5], displays an analogue behavior once computed for our 
optical lattice scheme. 

To compute p (E), we consider the Fourier transform of the 
free Hamiltonian Eq. 2 with periodic boundary conditions 

a,o'=l J 7=1 

where C aa ' = 8°- {a+l) + 5 <J ' (cT+1) (the N+l layer is iden- 
tified with the first one due to the periodic boundary condi- 
tions) is the the matrix that generates the cyclic group Z#. 
This Hamiltonian can be diagonalized in the space of layers, 



dx 


5 Er 


10 Er 


15 Er 


20 Er 


40 Er 


0.05a 


. 93625 


. 903827 


.879944 


.860432 


.803103 


0.1a 


0.768873 


. 66791 


. 60015 


.548733 


.416638 


0.15 a 


0.554949 


.404605 


.318226 


.260268 


. 140275 


0.2a 


. 353203 


.2018 


. 131853 


.0923365 


0.0308881 


0.25 a 


.19915 


.0833781 


. 0429945 


. 0246996 


0.00449657 



Table I. The ratio Ru between the inter-layer and the onsite interac- 
tion terms for D = 3, as a function of the optical potential V and of the 
x-component of the displacement d x = |d|/\/3- As confirmed by the 
numerical results, Ru is monotonically decreasing function of V and 
d x . Only for V = 5Er, Ru is never below 1/10 in the range of d x com- 
patible with a bivolume hopping term, see figure 2. The inequality 
Ru < 1/10 for d x /\/3 > a/4 and V around IOEr and higher, ensures 
that a four 3-dimensional layer with negligible inter-layer interaction 
can be consistently achieved for such values of the potential. For V 
greater that 20£r a five 3-dimensional layer is possible. 
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H=Z a=hN J^EW(k,e)a^aP with 

E D > N (k,o) = -2J £ cos&;-2/'cos(^% 

j=l,D 

where the tower of Kaluza-Klein modes modify the disper- 
sion relation introducing terms proportional to the coupling 
between layers J'. It follows that the density of states is 




First, we observe that for N 1, the sum can be approximated 
by the integral over kr> + \ = 2n(<3 — 1)/N, and Po,n(E) — > 
Pd+i,\{E). Second, Po,n(E) can be computed analytically in 
the low-energy limit, E ~ E m j„ = —2(DJ+J'), for / = /', or 
for/' <C J, by taking the continuous limits cosfcy — s- (1 — jlcj), 
j= 1,...,D + 1, or j = 1,...,D, respectively. In the latter case, 
the density of states for N = 2, in terms of the one of two 
uncoupled layers i.e. Pd.n=\ x j( ~ J mh ) a / 2 ~ 1 , is 

Pd.n=2~ 3Pai((1 + E=^) T_1 + (1-f4-) 1_1 ) 

= Pd,i(1 + ^ — — ( £ _ £ . j +•••)• 

In particular, the above expression provides an interesting ex- 
perimental signature for a bivolume, D = 3. 

A many body observable: MI-SF transition. In the pres- 
ence of interactions, the effective dimensionality of the sys- 
tem can be experimentally detected by measuring the location 
of the Mott-insulator-to-superfluid phase transition. For sim- 
plicity, we focus on the bosonic case. It is well known that 
the ground-state of Bose-Hubbard Hamiltonian in any dimen- 
sions presents two phases: for J/U -C (J/U) c , the ground- 
state is a Mott-insulator with definite local occupation. For 
J/U 3> (J/U) c the ground state is a superfluid state with 
all the atoms condensed in the single-particle state with null 
pseudo-momentum. The two phases, Mott-insulator and su- 
perfluid, are characterized by the local order parameter (a q ). 
As the Mott-insulator ground-state has a well-defined local 
occupation we have (a q ) = everywhere, i.e. no atom num- 
ber fluctuations. For the superfluid state (a q ) ^ 0, as the 
occupation is local in momentum and not in position. The 
critical point (J /U) c separates the phases with zero and non- 
zero order parameter, or equivalently, where the symmetry 
fl q — > e'^flq is broken and where it is not. This critical point 
(J /U) c depends on the dimensionality of the lattice. 

To be more precise, the effective Hamiltonian of a multi- 
state optical lattice in D spatial dimensions corresponds to 

H = -£(f>£if.4 CT) +J , a i T a+ ^al a) ) +H.C.+ 
r,a j=\ 

+ ^ni a) ((n { r a) 1) +Ru(n { ° +1) - 1) -2£) , (7) 

r,cr 

where U = i^L Jd D r \w(r)\\ n [ r a) = o^V 5 is the num- 
ber operator, and £i = jj is the chemical potential measured 
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Figure 3. Ground-state order parameter \jf as a function of J/U for 
different values of coupling j' between the two volumes. For J 1 = 
we obtain the same result as for D = 3. As the coupling between the 
layers increases, the critical point (J/U ) c approaches its D = 4 value 
for J' = 2J, which corresponds to j' = J, when periodic boundary 
conditions are chosen. The inset shows the critical point (J/U) c as a 
function of the interlayer coupling. 

in units of U. If Rjj <C 1 the above Hamiltonian is the Bose- 
Hubbard Hamiltonian in D + 1 dimensions. The value of Rjj 
is controlled by the the depth of the optical potential V, and 
the inequality is always satisfied for large enough V. Thus, we 
can appropriately tune the scattering length a such that MI-SF 
occurs for a very large value of V. 

In Fig. 3 we plot the order parameter computed using a 
Gutzwiller ansatz as a function of J/U for a lattice made of 
two D = 3 layers. Different curves correspond to different 
values of the coupling between layers, As / — > 2J, which 
corresponds to J 1 — > J in the case of periodic boundary con- 
dition, the critical value approaches the known value for four 
dimensions, in mean-field theory. The inset shows the theoret- 
ical prediction of the shift in (J/U) c as the coupling between 
the two layers increases. 

Within the mean-field approximation, the value of (J/U) c 
can be computed analytically extending the second order per- 
turbative approach of [34, 35] to the Hamiltonian in Eq. 7 for 
N layers periodically identified and Rjj = 0. Indeed, due to the 
periodic boundary conditions in the compact extra-dimension 
-the lattice is assumed to be sufficiently extended in the other 
D dimensions such that boundary conditions do not matter- 
the order parameter is constant (and it can be taken real) on 

the lattice \j/ = (flr ), Vr, rj. The critical value is found when 
the symmetric phase y/ = becomes unstable, i.e. when the 

- IHjBI — o. Such quantity can be computed exactly, 

within the mean-field approximation, by treating the hopping 
term as a perturbation at second order. The result is 



aU c = 2n+ \J{2n+ 1) 2 -1, (8) 
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where a = D/(D + J'/J), U c = U/{2DJ) is the critical value 
made independent of the connectivity of the hypercubic lat- 
tice, and h is the occupation of the Mott state in each species. 
For /' = 0, it reduces to the known expression. Let us stress 
that Eq. 8 is not affected by the number of species, N. In 
fact, for/' = / it follows that (J/U) c for D-dimensional model 
with N layers coincides with critical value for the D + 1 -Bose- 
Hubbard. This the case because in the mean-field approxima- 
tion only local properties like the number of neighbors count, 
and nothing can be said about global properties like topology 
and size of one spacetime direction. Although this prediction 
cannot be exact -ra-point correlation functions for m > N cer- 
tainly distinguish whether a spacetime direction is compact- 
ified on a circle or not- we expect that the mean-field result 
is qualitatively correct, and that a small number 3D layers are 
sufficient to reproduce the AD critical behavior. 

Conclusions. We have presented a strategy to produce a 
quantum simulation of an extra dimension. The recipes we 
have proposed to engineer 4D models pave the way to the 
study of novel phases, not accessible in 3D condensed mat- 
ter world. Furthermore, together with the capacity to simulate 
the propagation of pseudoparticles in nontrivial gauge fields 
and in a gravitational background [36], the present proposal is 
a step to a complete tool-box for simulating quantum field the- 
ory scenarios, in and beyond the Standard Model of particle 
physics. 
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